## # A tibble: 277 × 4
##    folder           plot_id       simul         simulation_id
##    <chr>            <chr>         <chr>                 <int>
##  1 output-cmd_1.txt RETZ_00972_01 simulation_10             1
##  2 output-cmd_1.txt RETZ_00972_01 simulation_11             2
##  3 output-cmd_1.txt RETZ_00972_01 simulation_12             3
##  4 output-cmd_1.txt RETZ_00972_01 simulation_13             4
##  5 output-cmd_1.txt RETZ_00972_01 simulation_14             5
##  6 output-cmd_1.txt RETZ_00972_01 simulation_15             6
##  7 output-cmd_1.txt RETZ_00972_01 simulation_16             7
##  8 output-cmd_1.txt RETZ_00972_01 simulation_17             8
##  9 output-cmd_1.txt RETZ_00972_01 simulation_18             9
## 10 output-cmd_1.txt RETZ_00972_01 simulation_19            10
## # ℹ 267 more rows
## # A tibble: 278 × 4
##    folder           plot_id       simul         simulation_id
##    <chr>            <chr>         <chr>                 <int>
##  1 output-cmd_1.txt RETZ_00972_01 simulation_10             1
##  2 output-cmd_1.txt RETZ_00972_01 simulation_11             2
##  3 output-cmd_1.txt RETZ_00972_01 simulation_12             3
##  4 output-cmd_1.txt RETZ_00972_01 simulation_13             4
##  5 output-cmd_1.txt RETZ_00972_01 simulation_14             5
##  6 output-cmd_1.txt RETZ_00972_01 simulation_15             6
##  7 output-cmd_1.txt RETZ_00972_01 simulation_16             7
##  8 output-cmd_1.txt RETZ_00972_01 simulation_17             8
##  9 output-cmd_1.txt RETZ_00972_01 simulation_18             9
## 10 output-cmd_1.txt RETZ_00972_01 simulation_19            10
## # ℹ 268 more rows

complete_scene <- import_output_scene( output_name = “complete”, base_path = “C:/Capsis4/data/forceps/clementine/Test_protocole/”, output_file = “output/protocole_completeScene.RData”) %>% filter(date <= “2100”) %>% group_by(plot_id, date, scenario) %>% summarise( totalVolumem3.ha = sum(totalVolumem3.ha, na.rm = TRUE), totalBiomass_t_ha = sum(totalBiomass_t_ha, na.rm = TRUE), .groups = “drop” )

Etats initiaux des placettes

État initial des placettes
Plot ID Surface Type de peuplement Essence déterminante Essence dominante 1 Essence dominante 2 Structure et occupation du sol Surface terrière (m2/ha) Diamètre moyen (cm) Âge médian (ans)
RETZ_00230_03 2.50 FHETE HET HET F 0.0 7.5 12.0
RETZ_00816_02 0.94 FHET1 HET HET F 17.5 12.5 12.0
RETZ_00901_01 10.23 FHETG HET HET CHP F 21.0 57.5 144.5
RETZ_00972_01 4.59 FHETG HET HET CHS F 31.0 57.5 144.5
RETZ_01008_05 3.38 FHETS HET HET CHS F 0.0 3.5 12.0
RETZ_01067_02 1.73 FCHPE CHP CHP F 0.0 7.5 12.0
RETZ_01210_04 3.10 FCHSE CHS CHS F 0.0 7.5 12.0
RETZ_01307_02 1.28 FDOUG- DOU DOU F 29.0 52.5 59.5
RETZ_01368_01 7.91 FHETM HET HET CHA F 31.0 37.5 114.5
RETZ_01415_02 2.23 FCHPG+ CHP CHP CHA F 17.5 62.5 NA

Exploration des itinéraires (pour validation des dynamiques forestières)

Pour rappel des itinéraires :

Je ne suis pas sûre de pourquoi il y a des chutes de surface terrière dans la simulaion 3 (sans coupe), l’itinéraire étant bien : 150_3_0.5_0%_FSyl-80. peut-être juste mortalité naturelle

Calcul et visualisation des indicateurs par placette

Sur les 80 ans de simulations on calcul les indicateurs suivants :

Indicateurs globaux

Relations entre indicateurs globaux

A faire